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Abstract. 

We present a brief review of the polarization properties of the cosmic mi- 
crowave background in dark matter models for structure formation. Quite 
independently of the model parameters, the polarization level is expected 
to be ~ 10% of the anisotropy signal at angular scales < 1°. Detections 
of polarization at larger angular scales would provide a strong evidence in 
favour of an early reionization of the intergalactic medium. 



1. Introduction: some historical remarks 

Most of the early theoretical work on the polarization of the Cosmic Mi- 
crowave Background (CMB) was focused, after Rees pioneering work [1], 
on anisotropic cosmological models [2,3,4,5,6,7,8,9]. The degree of Faraday 
rotation expected in the presence of an universal magnetic field and the 
use of polarization measurements to constraint the amplitude of such field 
were also considered [14,33]. More recently, it has been shown that even 
in isotropic cosmological models the anisotropic component of the CMB is 
polarized [10]. Detailed numerical predictions have been made for dark mat- 
ter dominated models with adiabatic fluctuations [see e.g. 11,12,13], with 
and without an early reionization of the intergalactic medium [15]. These 
calculations have shown that the level of polarization can be 10 percent 
of the anisotropy signal. After the COBE/DMR result [16], new attention 
has been dedicated to the tensor modes of metric fluctuations, which also 
produces anisotropy on large angular scales [17,18,19,20,21,22,23]. The po- 
larization due to a background of primordial gravitational waves has been 
widely discussed [24,25,26,27,28]. For describing the statistic of the polar- 
ization field was also introduced the polarization - anisotropy correlation 
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Figure 1. Thomson scattering of a photon by an electron. 



function [29,30], while other authors [31,32] have shown that neglecting 
polarization yields a theoretical overestimate of the anisotropy at small 
angular scale. 

From the experimental side, in spite of a continuous increment in the 
experimental sensitivity, no polarization was found and only upper limits 
were given [34,35,36,37], with the best upper limit to date of ~ 25/iK from 
the Saskatoon experiment [38]. As we show in Section 5, the level of CMB 
polarization expected is in most of the models at least a factor 10 below this 
limit, so is not clear if the present sensitivity of the CMB experiments is 
sufficient to detect polarization. However, in view of forthcoming high sen- 
sitivity new experiments, it is of interest to discuss the general properties 
of the polarization pattern and its dependence to the various cosmological 
parameters. So, the aim of this work is to review the basics steps behind 
the theoretical calculation of CMB polarization. The plan of the paper is 
as follows. In Section 2 we briefly review the definition of the Stokes pa- 
rameters and their variations in a Thomson scattering. In Section 3 and 4 
we write the set of equations necessary to describe anisotropy and polar- 
ization of the CMB. In Section 5 we review some of the results obtained by 
numerically integrate those equations. Finally, in Section 6, we summarize 
the main findings. 

2. An Elementary Description of the Polarization of Light 

For an elliptically polarized wave, the components of the electric field along 
two orthogonal directions, £ and r say, can be written as : 
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E{ = E^sm(ujt - ei) 
E T = sm(ujt - e 2 ) 



(1) 



where E® T and are constants. The polarization of the radiation field is 
conveniently described in terms of the Stokes parameters : 
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The parameter / is proportional to the intensity of the wave (we omit 
the proportionality factor) V measures the ratio of the principal axes of 
the polarization ellipse while Q or U measures the orientation of the ellipse 
relative to the £ axes. In general I 2 > Q 2 + U 2 + V 2 , the equality holding 
for an elliptically polarized wave. 



A clockwise rotation by an angle S of the £ — r axes in the polarization 
plane leaves unchanged the I and V parameters and it is equivalent to 
apply the operator: 
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(3) 



to the vector / = (i^, I T , U, V) . 

In the Thomson scattering, the light scattered in a direction making an 
angle with the direction of incidence (see Figure 1) is 



E l = \Jl aTE l cos sin M - e i) 
E r = \Jl^rE Q T sm{ut- e 2 ) 



(4) 



In analogy with the equation (2) the Stokes parameters of the scattered 
light are: 



J| = (3/2)a T / 5 cos 2 e 
J* = (3/2)a T I T 
U s = (3/2)CT T [/cose 
I V s = (3/2)<7 T Fcose 



(5) 



or, in matrix form, I s = axR x /, where 




Figure 2. Coordinate system needed to describe Thomson scattering in the Lab frame. 
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In order to study the variations of Stokes parameters in the Lab frame 
(see Figure 2) we have to : 

1. apply the transformation L(—i\) to /, where i\ is the angle between the 
meridian and the scattering planes. In this way we obtain the Stokes 
parameters of the incident light in the frame of Figure 1. 

2. apply R to these parameters in order to obtain the Stokes parameters 
of the scattered light, again in the frame of Figure 1. 

3. apply the transformation L(tt — {2), where i% represents the angle be- 
tween the plane OP2Z and OP\P2- In this way we are back to the Lab 
frame. 

Thus, the radiation scattered in the (0, </>) direction, relative to the Lab 
frame, can be written as [39]: 

I s (9, 4>) = -L / [P(6, cf>; 6', </>') x 1(6', <f/)]dtf (7) 

47T Jilt 
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where 



P = Qx [P°(ji, //) + J (I - ^)J(1 - v' 2 )P\^ 4>, A 4 ') 4>') + P 2 (^ <t>, M 
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and /i and // are denned as cos6> and cos#', respectively. 
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3. The Boltzmann Transfer Equation for Polarized Light 

In order to study anisotropy and polarization of the CMB we need to write 
down the transfer equation for the Stokes parameters. We restrict ourselves 
to isotropic universes where, to zero-th order, anisotropy and polarization 
vanish. The perturbations to the Stokes parameters and to the other rele- 
vant quantities (see section 3.3) are written in the synchronous gauge for- 
malism. Following Peebles [40,41] we introduce the fractional fluctuations 
of the Stokes parameters as follows: 
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where (i,q,u,v) are functions of the observer position x, of the line of sight 
observation 7 = (71,72,73) and of the cosmic time t. To first order the 
transfer equations becomes 
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where y = —2h a p~f a jp is the term containing the linear perturbation to the 
metric tensor, (l s ,q s ,u s ,v s ) are evaluated in the comoving frame and refer 
to the radiation scattered in the 7 direction, and a is the scale factor. In 
order to avoid spatial dependence it is convenient to work in Fourier space. 
We choose for each k mode a reference system with the z axis parallel to 
k, in order to achieve an azimuthal symmetry. 

3.1. SCALAR MODES 

For scalar modes, the only non vanishing components of the perturbed 
metric tensor are the diagonal ones : hn, ^22, ^33 (^00 = because of the 
chosen gauge). Thus, the gravitational term in equation (14) has the form: 
y = (1 — 3fi 2 )Ji33 — (1 — fj> 2 )h, and each Fourier mode is independent of 
the azimuthal angle 4>. After integrating over this angle equation (7) it can 
be proved that P 1 and P 2 give no contribution. Therefore we can assume 
U = in this case. Also the equation for V is decoupled from the others: 
if V vanishes at the beginning, it also vanishes afterwards. Therefore only 
the perturbations 1 and q of the / and Q parameters are of interest. Their 
evolution is described by the following transfer equation [10,11,12,13]: 

Ko^a)-©— a>^(^*'-c + r)) 

(15) 

where y = (l — 3fj, 2 )h^ — (l — fj 2 )h is the term taking into account the effects 
of gravitational potential, and where the 2x2 matrix M$ is composed by 
the first two rows and columns of the matrix P° in the (I, Q, U, V) basis: 

3-^- M 2 + 3^y l-/i' 2 -3/i 2 (l-/i' 2 ) \ 
1 - 3/Z 2 - fi 2 + 3^V 2 3 - 3,u /2 - 3/i 2 (l - fi' 2 ) ) [lb) 

The angular dependence in equation (15) can be eliminated by expand- 
ing 1 and q in Legendre polynomials : 
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* = E(4W P 2^) + ial e+1 (t)P 2e+1 (fi)) (17) 

Because of the orthogonality of the Legendre polynomials, equation (15) 
becomes: 

1 + ^-2/5 = <rrn e (ao - 4/w 6 - i + P 2 (/x)l/2(a 2 /5 - % + r? 2 /5) 

% + i <fq = a T n e (-q + 1/2(1 - P 2 ( / x))(- < r 2 /5 + % - r? 2 /5) 

(19) 

These equations are coupled together through the quadrupole term, i.e. 
the radiation must have a quadrupole anisotropy to get polarized. 

3.2. TENSOR MODES 

For tensor perturbations the only non vanishing components of the per- 
turbed metric tensor are h\\ = /i 22 = h + and hyi = /i 2 i = h x where the 
two values h + , h x refer to the two polarization states of the gravitational 
waves. The equation of transfer has the following form [24] : 




(20) 

where now y = —h + (l — y?) cos(2</>) + h x (1 — y?) sin(2</>), and the 3x3 
matrix Mt is composed by the first three rows and columns of the matrix 
P 2 in the (J, Q, U, V) basis: 



3 /K_(/i)if_(//)cosA^ -K_{n)K + (n')cos^ -2/i / X_(//) sin 
Mt = — I K + (n)K-(n ) cos -K + (p)K + (ji') cos A^ 2fi'K + (p) sin A 

\ fj,K-(n') sin A^ —fiK + (p')s'mA^, 2////cosA (? j / 

(21) 

with K±(fi) = 1 ± fi 2 , A^ = 2(<^>' — <^>). The particular form of the metric 
tensor makes i still dependent on <p in spite of the choice of a special refer- 
ence system. However, this dependence is not too cumbersome. In fact, it is 
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possible to introduce a change of variables [24]_ to eliminate the dependence 
on the azimuthal angle. The new quantities, /, Q and U are related to the 
old ones by the following relation: 

(j)) = - n 2 ) cos 20 + I x - p 2 ) sin 2<j) 

Q((j.,<j>) =Q+(^)(l + / t/ 2 )cos2</» + Q x ( / u)(l + /[ x 2 )sin20 (22) 
U(p, 4>) = -U + 2p sin 2<j) + U x 2p cos 2<j) 

It is easy to prove that, with these new variables, only P 2 provides a non 
vanishing contribution to the integral of equation (7) over (p. This is why 
we have considered only this term in equation (21). Also, as the Boltzmann 
equation becomes independent of (ft, we can still develop fluctuations of I,Q 
and U in Legendre polynomials. Thus, equation (20) becomes: 

' 4 + iJ fk ~ 2h e = -a T n e (l e + 

< 4 + i^q t = -a T n e (q e - *) (23) 
, q e + u e = 

where 



* = 3/5^,o + 6/35^,2 + l/2107? e ,4 - l/10er e , + l/35a e , 2 - l/2Wa e , a (24) 

and e identifies either the + or the x polarization state of the gravitational 
wave. 

3.3. NUMERICAL CALCULATIONS 

We restrict ourselves to a Universe composed by baryons, cold dark matter, 
photons and three families of massless neutrinos. The equations describing 
anisotropy and polarization of the CMB have been written above. In Fourier 
space, the equations describing fractional fluctuations in the remaining cos- 
mic components are [42,43,44]: 
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v b + H(t)v b = a T n e /, 
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for scalar and tensor fluctuations, respectively. 

Eq.(25) describes the fluctuations in the massless neutrinos. We follow 
this component only when the perturbation proper wavelength is larger 
than one tenth of the horizon. Afterwards, free streaming rapidly damps 
fluctuations in this hot component. 

The time evolution of the baryon and CDM density contrasts and of 
the baryon peculiar velocity are described by Eq.(28), (30) and (29) respec- 
tively. The system for the scalar fluctuations is closed by Eq.(26) and (27) 
describing the field equations for the trace and the 3 — 3 component of the 
metric perturbation tensor, while Eq.(31) is all we need to describe the evo- 
lution of the metric perturbations for tensor fluctuations. We numerically 
integrate the previous equations from redshift z = 10 7 up to the present. 
The abundance of free electrons, n e , is evaluated following a standard re- 
combination scheme [45,46] for H and 4 He, taken in the ratio 77 : 33. In 
the following we also consider the possibility that the universe reionized in- 
stantaneously at redshift z r h « 1000, and remained completely reionized 
up to the present. 



4. Computing the Correlation Function for Anisotropy and Po- 
larization 

Under the assumption of gaussian initial conditions, the statistical proper- 
ties of the CMB anisotropy and polarization patterns are fully described in 
terms of their correlation functions. The stochastic anisotropic component 
of the CMB is conveniently expanded in spherical harmonics: 5T{^)/T$ = 
J2im a lmXm{l) ■ The coefficients a\ m are random gaussian variables with 
zero mean and rotationally invariant variances, Ce = (\ a\ m | 2 ). The mean 
(over the ensemble) correlation function of the anisotropy pattern has the 
standard expression: 




(32) 




Figure 3. Laboratory reference system. 



where cos 9 = 71 • 72 , and 



Here the primordial power spectrum of scalar fluctuations is assumed 
to have the standard form P(k) = A$k ns . In the case of tensor fluctua- 
tions, the change of variables needed to achieve rotationally symmetry [see 
equation (22)] must be taken into account. Then, the correlation function 
of the CMB anisotropy induced by tensor modes reads: 



,ar(7i)£T(7 2 ). Ai 



N T T ' 128vr 3 
where 



(34) 



= K -{^' i)K-{^' 2)^11^2* cos(20i ) cos(2</> 2 )+cr£cr^* sin^i ) sin(20 2 )] 

(35) 



where the power spectrum of metric fluctuations due to tensor modes is 
assumed to be P(k) = j4r/c nT ~ 3 . Assuming = <5^ x and making some 
algebraic manipulations yield: 
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(36) 

where 

T = [2(7x • 72 - AM? ~ K^\)K^' 2 )] (37) 
and, finally [47,48,49], 

JO ^0 



with 



and 



c, - A ^ e+2 y- r '^ (fc) iv -^ (39) 



nj (2^-l)(2£-3) (2£-l)(2* + 3) (2£ + 5)(2£ + 3) 1 ; 

As discussed in Section 2, the Q and U Stokes parameters vary because 
of a clockwise rotation S of the reference system in the polarization plane : 

Qs = Q cos(2~ ) + 17 sin(2H) , . 

U s = -Q sin(2 ~ ) + U cos(2~ ) 1 j 

So, the perturbations to Q and U in the Lab frame are related with 
those in k space by the following relation: 



Q(x,9,4>) 1 



T 32vr 3 
U(x,0,</>) -1 



J g(fc)e ife cos[2~(fc)]d 3 £; (42) 
J g(fc)e ife sm[2~(A?)]d 3 fc (43) 



T 32vr 3 

The correlation function for the Q parameter for scalar modes can be writ- 
ten as follows: 



( ^^TT ) = ^ VeiVto cos(2Hi) cos(2E 2 )Pn(A 1 )Pe2(A2)k ns d 3 k 

(44) 
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The correlation function for U has a similar expression with cos2H — > 
sin2H. Let us identify the line of sight 71 with the z-axis of the Lab frame. 
In this case, (#i,</>i) = (0,0), 6^ = 9/; and <^ = — ^ (see Figure 3). In the 
small angle approximation, cos(2Hi) ~ cos(2H2) ~ cos(2<^), and equation 
(44) yields : 

m9m )=m+B(eA) (45) 

where 

m = 25^3 /EvUk)rimPn^ k )Pe2(^2)k ns d 3 k (46) 

and 

5(M) = ^3 / E^(fcK 1 (A;)P, 1 ( /U ' fc )P, 2 ( / u' 2 )cos(4^)^ Q! 3 A; (47) 
The first term has the standard expression : 

A W = 4^ E( 2 ^ + l)Q Q ^(cos 9) (48) 



with 



' lOvrio (2^ + l) 2 1 j 

For the second term, we can develop Pa (cos 0^) i n associated Legendre 
polynomials with m = 4: 

00 

Pn(cos(0 fe )) = ]T a/ M 'i^(cos((9 fc )) (50) 

where 

2f + l(f-4)< 
a/W = 2 (f + 4)T A <M' ( 51 ) 

and = P^Pfid^osOk) has the following values : 
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A e>i > =8(£' 2 +£' -3(£ 2 + £ + 2)) £> = £ + 2,£ + A, ...,£ + 2n ^ 

Ae,e = £' = £+ 1,£ + 3, ...,£ + 2n + l 

=0 f < £ 

From the definition of the spherical harmonics it also follows: 
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Pf,{cose k ) ■ cos(40 fc ) = ^ 2/+! (^4)] + 4 $) ( 53 ) 

Now, execute the following steps: 

1. Insert the (50) in the (47) 

2. Develop the product between the Legendre polynomial and the cosine 
with the (53) 

3. Integrate in d£lk 

4. Use again the (53) to transform the spherical harmonics in Pf 

5. Develop the Pf in Pj>. 

At the end, in the small angle approximation, it is possible to write [29,50,71]: 



to ^0 



where the two terms and Be are defined, for scalar perturbations, as 
follows: 

C ? = 7FZ ( 7^* nS+2 dk (55) 



16tt J (2^ + l) 2 

Be = £g J^ A '* A ™ J *V s4 '* < 56 > 

For tensor fluctuations, the calculation is similar. The final result is 
[30,50,71]: 

{9 n l9 n 1) = ^E^ + ^c? + cos^)B e )p £ (co S e) (57) 
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(2£-l)(2£-3) 
(61) 



and ^1^2 = jPa(^)^2 4 (^)^- 

An interesting case is the correlation function between CMB anisotropy 
and polarization. For scalar fluctuations and in the small angle approxima- 
tion, the result is [29,50,71]: 
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5. Numerical Results and Discussion 

With the formalism developed in the previous Sections, we are now able 
to make theoretical predictions for CMB anisotropy and polarization. As 
stated before, we restrict ourselves to the Cold Dark Matter scenario. This 
is not quite enough to completely define the model, as we have to deal with 
quite a number of parameters. We have to fix : i) the total density pa- 
rameter, Qq, and the cosmological constant, A; ii) the baryonic abundance, 
f2&; iii) the Hubble constant; iv) the primordial spectral index for spectral 
fluctuations; v) the relative amplitude of scalar and tensor fluctuations; vi) 
the spectral index for tensor fluctuations; vii) the thermal history of the 
universe; viii) the overall amplitude for scalar fluctuations. 

The baryonic abundance is quite severly restricted by primordial nucle- 
osynthesis. We consider the fiducial value of fifc = 0.05 ±0.02 as representa- 
tive of the possible uncertainty in this parameter. We remind that changes 
in fife yield variations in the pressure of the photon-baryon fluid before re- 
combination, and then variations in the amplitude of the first acoustic peak 
in the anisotropy power spectrum. 
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For flat models the value of the Hubble constant is usually taken to be 
Ho = 50Kms~ l /Mpc for age considerations, even if the estimated globular 
cluster age allows for slightly different values: AOKms' 1 /Mpc < Hq < 
QhKms' 1 /Mpc [63]. Small variations in the Hubble constant yield huge 
variations in the radiation power spectrum, modifying both the amplitudes 
and positions of the acoustic peaks. 

The primordial index for scalar fluctuations, ns, is usually taken to be 
unity, as inflation suggests. However, in more general inflationary scenarios, 
ns can be either smaller or larger than unity [64]. This modifies the relative 
power between the anisotropy on small and large angular scale. Power-law 
inflationary models predict ns < 1, but also a background of gravitational 
waves. A standard prediction is that the ratio of the quadrupoles induced 
by scalar and tensor fluctuations is: 

^ = -7n T ~7(l-n s ), (65) 

which allows to relate amplitudes and shapes of primordial power spectra 
for scalar and tensor fluctuations, respectively [18,49,64]. 

The thermal history of the universe, in its standard form, assumes re- 
combination of the primordial plasma at redshift ~ 1000. However, both 
the Gunn-Peterson test [65] and the enriched composition of the intraclus- 
ter medium [74,75] suggest the possibility of a considerable energy release 
during the early stages of galaxy formation and evolution. 

Finally, the amplitude of fluctuations it is still unknown from first prin- 
ciples, and it is fixed in order to match the observed rms temperature 
fluctuations ( 29 ± luK) of the COBE/DMR anisotropy maps [53]. 

In Figure 4 we show theoretical predictions for CMB anisotropy and 
polarization of a standard Cold Dark Matter Model with Qq = 1, = 0.05, 
Hq = BOKms^ 1 j Mpc, ns = 1 and standard recombination. The anisotropy 
power spectrum for scalar fluctuations has a flat behavior at low fs, where 
the Sachs- Wolfe effect [66] dominates, and a structure of peaks at higher €s, 
defined by the acoustic oscillations in the photon-baryon fluid experienced 
before recombination by fluctuations smaller then the acoustic horizon. 
The damping at high £'s is due to the finite thickness of the last-scattering 
surface [67]. The first peak at £ ~ 200 corresponds to fluctuations that 
entered the horizon at recombination, which subtends an angle ~ 2°h~ 1 . 
The polarization power spectrum has instead power only at £ > 200, i.e. 
on scales < 2°h~ l . The case of pure tensor fluctuations is show in Figure 
4 only for the didactic purposes. In this case, the anisotropy spectrum has 
power only at £ < 200 and the polarization spectrum shows a prominent 
peak at I ~ 100. Note that the polarization spectrum has an amplitude 
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Figure 4- Anisotropy (top) and polarization (bottom) power spectrum for scalar (left) 
and tensor (right) fluctuations. Each model is normalized to COBE er(10°) = 



lower than the anisotropy spectrum by a factor ~ 10 2 and ~ 10 for scalar 
and tensor fluctuations, respectively. 

Real experiments are sensitive to a limited region of the power spectrum, 
because of the antenna beam and modulation techniques. For anisotropy 
experiments, this effect can be parameterized by a window function, We, 
so that the variance of temperature fluctuations detected by an experiment 
can be written as: 



(66) 



A similar expression holds for the variance of the fluctuations of the Q 
parameter: 
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(67) 
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Figure 5. R.m.s. values for anisotropies (top) and polarization (bottom) scheduled 
for scalar (left) and tensor (right) models for different experiments [53... 61]. The tensor 
model has nr = and the scalar model has ns = 1. Each model is normalized to COBE 
<r(10°) = 29fiK. The dot of the left side represents the values for a scalar reionized model 
at z r ~ 70. 



In fact, it can be proved that the azimuthal contribution [see equations 
(54) and (57)] to the Q variance vanishes. In equation (67) we use the 
same anisotropy window functions, in order to give an order of magnitude 
estimate of the level of measurable polarization at different angular scales. 

In Figure 5 we plot the expected rms values for CMB anisotropy and 
polarization using 15 different window functions corresponding to 9 differ- 
ent anisotropy experiments. The level of polarization from scalar modes 
is below the current experimental sensitivity, even for small scale experi- 
ments such has Saskatoon or CAT that are sensitives to multipole £ ~ 400 
where the polarization has the first two peaks. The rms level from pure 
tensor modes, even for the MAX experiment that seems to have the best 
window function, is below 0.5 fiK, so the separation between scalar and 
tensor fluctuations do not seems to be at hand with polarization measure- 
ments [25,26,27,28,73]. This can be done by combining anisotropy mea- 
surements at both large (where tensor modes could contribute) and small 
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Figure 6. Dependence of polarization on cosmological parameters. 

(where tensor modes do not contribute) angular scales [49,72,73]. An accu- 
rate mapping of the anisotropy pattern with both high sensitivity and high 
angular resolution will be provided by planned, dedicated space missions 
such as COBRAS/SAMBA [68] and MAP [52]. At the moment the bulk of 
degree-scale detections, combined with the COBE/DMR and Tenerife ex- 
periments, seems to suggest a spectral index for scalar fluctuations n$ > 1 
[69] and a negligible contribution of tensor modes. 

As mentioned above, there are several free parameters, each with its 
own uncertainty, which define a theoretical model. So, it is interesting to 
explore the sensitivity of the polarization level relative to the anisotropy 
one. To show this, we plot in Figure 6 the quantity /Q as a function 
of £, for different choices of the model parameters. Generally speaking the 
effects of the variation of these parameters on anisotropy and polarization 
are similar: both quantities tend to decrease with increasing Hq and tend 
to increase when a cosmological constant Qa = 1 — is taken into ac- 
count. In particular, varying ns or fi& yields basically no variations on the 
Cf/Ce ratio. Moreover, decreasing the spectral index and adding gravita- 
tional waves increase the large scale polarization, but, as we have shown, 



19 



6,0x1 (T - 




(degrees) 



Figure 7. Polarization - Anisotropy correlation function. 



not enough to pass the threshold of present day detector sensitivity. So, 
even taking into account reasonable uncertainties in the parameters, its 
seems that only with an huge increment in the experimental sensitivity 
(see the accompanying paper by F. Melchiorri et al. in this volume) and/or 
a space mission [52,70,71] a robust detection of the polarization spectrum 
over a wide range of €s would be possible. Coulson et al. [29] suggested 
searching a correlation between the temperature in one direction and the 
polarization in a circle at distance G from that direction. The shape of the 
correlation function (62) measurable in this way is show in Figure 7. As we 
can see the cross correlation is positive on scales > 1°, negative on scales 
between 0.5° and 1° and positive again on scales < 0.5°. According to [52] 
the future MAP satellite would have the capability to measure the expected 
amplitude of this signal. 

The final item to be investigated is the dependence of CMB anisotropy 
and polarization on the thermal history of the universe. A reionization at 
z < 100 produces a new, later and thicker last scattering surface. The 
effect of such a new last scattering surface is to smooth the anisotropy on 
small angular scales and to leave unchanged the level of anisotropy on large 
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angular scales. The effect of reionization on polarization is to reduce the 
polarization on small scales but to increase the polarization level at large 
angular scales. This is shown again in Figure 6. Thus, possible detection 
of polarization between 1° and 10°, say, would be an evidence for an early 
re-heating of the intergalactic medium. 

6. Conclusions 

Numerical solutions of the Boltzmann transport equation show that a cer- 
tain degree of polarization must be present as a consequence of the primor- 
dial fluctuations responsible for the structure formation. The level of polar- 
ization depends strongly on the angular scale, much more than in the case 
of anisotropy, quite independently of the choice of the model parameters. 
At angular scales larger than one degree we do not expect detectable po- 
larization unless the universe was reionized at early times z > 40. At small 
angular scales the polarization may reach the 5 — 10% of the anisotropy. 
However, a polarization of a few percent at angular scales of 1°-10° could be 
explained only because reionization: a search for polarization at these scale 
is therefore important in the study of the thermal history of the universe. 
Also, it seems hard to disentangle tensor from scalar perturbations through 
measurements of polarization, due to the tenuity of the signal. 
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